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Abstract 

This paper revisits the problem of characterizing the highest spatial throughput - or spatial capacity - of 
wireless networks in [bits/s/Hz/m'^] in light of some recent results. We analyze the expected maximum achievable 
sum rates over a given area based on the capacity regions of Gaussian point-to-point codes in Poisson distributed 
wireless networks under two decoding rules, namely (i) IAN: treating interference as noise, and (ii) OPT: jointly 
detecting the strongest interfering signals, treating the others as noise. We prove that the average spatial capacity 
is the product of the network density and the expected value of the achievable coding rates that maximize the 
spatial throughput for different network realizations. Assuming that the closest interferer power approximates the 
aggregate interference caused by all transmitters that are treated as noise, we analytically derive several properties 
of the average spatial capacity and evince the better performance of OPT decoders. We also compare these results 
to the ones obtained for a scenario where transmitters code their messages at predetermined fixed rates that are 
tuned to optimize the average spatial throughput, regardless of particular realizations of the network. Using this 
approach, however, some rates may not be achievable, yielding outage events. We analytically show that, when the 
same decoding rule and network density are considered, the spatial-capacity-achieving scheme always outperforms 
the spatial throughput obtained with the best predetermined fixed rate strategy. We then discuss the validity of our 
approximation using simulation results and implementation issues related to different mobility patterns. 
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I. Introduction 

In the last few years, the demands for more efficient, reliable wireless systems induced network designers 
to think about alternative ways to supplement centralized cellular models. One interesting idea is to build 
a multi-layer network where macro-base-stations coexist with a great number of smaller cells, which in 
tum operate in a more distributed fashion (e.g. the concept of femto-cell networks [[I]). But differently 
from the centralized approach whose capacity are fairly well characterized by Shannon theory, the limits 
of distributed systems that work in interference-limited regimes are unknown except for few specific cases, 
as discussed in 

In 1978 Carleial formally stated the interference channel problem using arguments from information 
theory [|3l. Since then, several results and bounds have been proposed for the interference channel as 
discussed in f4, Ch. 6]. Although these works put some light on the problem, the capacity region of the 
simplest two- source-two-destination setting is still an open problem. Moreover, when multiples sources 
and destinations are considered, such capacity regions becomes even more elusive. 

Knowing such difficulties, some researchers have started investigating alternative approaches to better 
understand the limits of wireless networks with multiple communication pairs. Gupta and Kumar intro- 
duced in [|5|| the transport capacity metric to determine how many bits-meter a wireless network with 
uniformly distributed nodes can reliably sustain when its density grows to infinite (asymptotic analysis). 
After this milestone, many other papers have focused on a similar idea, finding the transport capacity 
scaling laws for different scenarios and under different assumptions. The monograph [6J compiles some 
of such studies. 

Franceschetti et al. presented another important result in |I71, where they applied an unconventional 
method to find the physical limit of wireless networks by using laws of electrodynamics. The authors 
further extended this approach in a recent work [IS] and determined the degrees of freedom of wireless 
networks based on electromagnetic theory. 

Nevertheless both Franceschetti's and Gupta's lines of research strongly rely on the asymptotic behavior 
of the network with the number of nodes infinitely growing, which may give an unclear picture of the 
physical or medium access control network layers' design. Bearing this aspect in mind, Weber et al. 
applied in BUl a statistical approach to characterize the throughput of wireless networks and then defined 
the transmission capacity as the highest spatial throughpu]^ achievable without exceeding a maximum link 
outage probability, using the density of active links as the optimization variable. An important aspect of this 

'in the literature, spatial throughput can be also referred to as area spectral efficiency 1101 . 
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work is the use of stochastic geometry (TT\ to characterize the node spatial distribution as a Poisson point 
process (PPP). Thereafter different strategies used in the wireless communications have been investigated 
such as interference cancellation, threshold transmissions, guard zones, bandwidth partitioning amongst 
others; the reference [[T2l compiles these results. 

In addition to them, we find in the literature other contributions using a similar approach. For example, 
Vaze studied in [|T3]| the throughput-delay-reliability trade-off in multi-hop networks using the metric 
random access transport capacity, which is an extension of the transmission capacity for multi-hop 
systems [[T2l Sec. 4.2]. In lfT4l the authors assessed how the hopping strategy affect the aggregate 
multi-hop information efficiency of ad hoc networks, while in [[TSl they studied a revisited version of 
the transmission capacity to compare different modulation -coding schemes. The work [[T6]| presents the 
transmission capacity optimization in term of the number of allowed retransmissions considering different 
MAC protocols. Ganti et al. generalized in flTTl the transmission capacity for different fading and node 
distributions for the high SIR regime. 

Apart from these recent papers that focus on the statistic quantification of the spatial throughput of 
wireless networks, the use of models from stochastic geometry dates back the early 80's, when Takagi and 
Kleinrock firstly introduced the idea of evaluating the aggregate interference power of Poisson distributed 
interfering nodes [[T8l . Thereafter, the subject has greatly developed and we can cite [fT9l - [|22l as relevant 
tutorials on how to apply stochastic geometry when analyzing wireless systems. 

This paper aims at characterizing the spatial throughput of wireless networks specifically when transmit- 
ter nodes are spatially distributed following a PPP. To reach this goal, we review the results due to Baccelli 
et al. [|23]| . where the authors studied the capacity regions of Gaussian point-to-point codes for interference 
networks. More specifically, we apply two different decoding rules described therein: treating interference 
as noise - the IAN rule - and joint detection of the strongest interferers ' messages and treating the others 
as noised - the OPT rule. Then, we use theorems from point process theory to obtain approximations and 
bounds of the highest achievable spatial throughput for a given network density, naming this maximum 
value as spatial capacity. 

For comparison purposes, we also evaluate the highest spatial throughputs achieved when all links 
use fixed (symmetric) coding rates (which is the most usual approach as in [[T2l . [|23l Sec. IV], [|24l '). 
We then proceed with an unconstrained optimization of the average spatial throughput over different 

^This rule splits the set of interferers into two mutually exclusive subsets: one contains the strongest interferers whose messages will be 
joint decoded with the desired one, and the other contains the transmitters with weaker detected power that will be treated as noise. This 
strategy is proved in 11231 the optimal for Gaussian point-to-point codes over interference channels, as discussed later on. 
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spatial realizations, where the optimization variable is the rate that the transmitters code their messages. 
Differently from the spatial-capacity-achieving scenario where the coding rates are tuned to be the highest 
achievable rates given the relative nodes' positions for each different spatial realization, the fixed rate 
scheme only cares about the average behavior of the network, resulting in decoding errors (outage events) 
for links whose capacity is below that predetermined rate. We analytically prove that, under the same 
assumptions, such a strategy performs worse than the spatial-capacity-achieving one. Our numerical results 
illustrate this difference as well as the advantages of using OPT instead of IA]v|^. 

We also simulate the network to validate our findings, justifying why our analysis is relevant even when 
our approximations are loose. Besides we discuss the feasibility of the decoding rules and optimization 
strategies for different mobility patterns: for (quasi-)static topologies, the spatial capacity can be achieved; 
for highly mobile topologies, it cannot and the fixed rate optimization with IAN turns out to be the most 
appropriate choice. From this perspective, we list the reasons to use cognitive, adaptive algorithms to 
improve the performance of wireless networks when mobility patterns are unknown. 

The rest of this paper is divided as follows. In Section HIl we first revisit the capacity region of 
Gaussian point-to-point codes over interference channels [|23l and then define the spatial capacity of 
wireless networks. Section III] presents the network modeling together with some general results for 
Poisson networks. Section |W] analyzes the IAN decoder, while the OPT is the focus of Section |Vl A 
comparison between the spatial-capacity-achieving and the symmetric rate scenarios is found in Section 
IVTl We discuss the tightness of our approximations and implementation issues in Section IVIIl followed 
by the final remarks in Section I Villi 

II. Capacity Region of Gaussian Point-to-Point Codes and Spatial Capacity 

This section reviews the capacity region of Gaussian point-to-point (G-ptp) codes for an arbitrary 
number of communication pairs, which was firstly stated by Baccelli et al. in ll23l Sec. II] and is also 
fundamental for the study carried out here. 

Let us assume a network with area A [m^] where K+1 source-destination pairs (also called transmitter- 
receiver pairs) coexist. Each source node i E [0, K] wants to transmit an independent message Mj G 
[l, 2"^*] to its respective destination i at rate i?, [bits/s/Hz], where n is the codeword length. Let Xj be 
the complex signal transmitted by source j E [0, K] and let Zi ~ CJ\f{0, 1) be the complex circularly 

'We do not assume any interference cancellation (IC) technique as in |12i Sec. 4.2], l25l -(27 1 since the OPT rule used in this paper 
always performs better than IC, as discussed in 0231 , ll24l . In addition, we consider no kind of strong coordination between the active links 
so that the interference alignment method 1281 is unfeasible; a detailed discussion about it can be found in 1231 Sec. III]. 
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symmetric Gaussian random variable that represents the noise effect at receiver i. The detected signal Yi 
at receiver i is then: 

K 

Yi = Y,9^JX, + Z,, (1) 

where gij are the complex channel gains between transmitter j (TXj) and receiver i (RXj). We assume 
that every transmitted signal is subject to the same power constrain of Q [W/Hz] so the received signal 
between TXj and RXj is constrained by Pij = \gij\'^Q. 

Each transmitter node uses a G-ptp code with a set of randomly and independently generated codewords 
x"(mj) = (xii, Xin){mi) following i.i.d. CJ\f{0,a^) sequences such that < cr^ < Q, where rrii E 
[l, 2"^'] , i E [0, K]. RXj receives a signal over the interference channel given by ([T} and then 
estimates the transmitted message as G [l, 2"^*] . An error event in the decoding happens whenever 

the transmitted message differs from the estimated one. Therefore the error probability of the G-ptp code 
is: 

1 ^ 

Next we use © to define the achievable rates and the capacity region for G-ptp codes. 

Definition 1 (achievable rates and capacity region): Let be the average of the error probability pn 
over G-ptp codes where n is the codeword length. Then, a rate tuple R = {Rq, Rk) is said to be 
achievable ifPn when n ^ oo. In addition, the capacity region using G-ptp codes is the closure of 
the set of achievable tuple rates R. 

This definition is important in our context because the spatial-capacity-achieving setting always requires 
achievable rates, as we will discuss later. But before that, we apply Definition [T] to establish the the capacity 
region of G-ptp codes. 

Theorem 1 (capacity region from /[23l/).- Let A be the set of all K + 1 transmitters in the network. Let 
Ai denote a subset of A that contains TXj with i E [0, K] and Ai its complement. RXj then observes a 
multiple access channel (MAC) whose capacity region TZi is computed as 



7^, = <j R: J^Rk < logs 
fee A 



y Ai^A}. (3) 



The capacity region TZ of the Gaussian interference channel with G-ptp codes is the intersection of the 
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capacity regions TZi of all TXj-RXj links with i E [0, K], i.e. 

K 

7^ = ^^- (4) 

i=0 

Proof: The proof of this theorem is found in [|23l Sec. II]. ■ 

The capacity region stated above requires a decoder that treats some of the interferers as noise, while 
others have their messages jointly decoded with the desired one. This result is the basis of the OPT 
strategy mentioned in the previous section and further studied in Section |Vl 

Now we define the network spatial throughput and then use Theorem [T] to define the spatial capacity 
as the highest achievable spatial throughput. 

Definition 2 (spatial throughput): Let A [m^] be the network area. Then the spatial throughput, denoted 
by S and measured in [bits/s/Hz/ni^], is defined as 

S = \Y.R^- (5) 

1=0 

Definition 3 (spatial capacity of a given spatial realization): The spatial capacity, denoted by S*, is 
defined as the maximum spatial throughput of the network such the rate tuple is achievable - R = 
{Ro,...,Rk) G Then, 

S* = max S. (6) 

Ren 

The spatial capacity reflects the highest sum of achievable rates over a given area and it may have a huge 
variation depending on the network topology (i.e. the number and/or position of the transmitter-receiver 
pairs). To deal with this issue, we opt for studying Poisson distributed networks that are analytically 
tractable, allowing us to derive approximate expressions for the expected spatial capacity S* (over different 
spatial realizations). 

III. Spatial Capacity of Poisson Networks 

Let $ be a two-dimensional homogeneous Poisson point process (PPP) with density A [nodes/m^] that 
characterizes the spatial distribution of transmitters (TXs) over M?. We assume that each TX is associated 
with one receiver (RX) located at a fixed distance d from it in a random orientationQ to establish a 
communication link; this is also known as Poisson bipolar model [[T9l . In addition, we consider that all 
TXs transmit information to their intended RXs over the same frequency band with width of W [Hz] 
(narrow-band) and using the G-ptp codes described in Section [III 

*Note that the RXs are not part of the process $. 
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For each realization of $, the network may have a different capacity region TZ and consequently 
a different spatial capacity S*. Even worse, when the network area is the infinite plane (i.e. R^), the 
capacity region given by equation dH) becomes impossible to be computed. Knowing these limitations, 
we define the average spatial capacity as the expected value of the highest achievable spatial throughputs, 
which allows us to analyze the performance of infinite networks over different spatial realizations of $. 

Definition 4 (average spatial capacity): Let R = {Rq, ...,Rk) be a tuple rate and IZ be the capacity 
region for a given network realization, then the average spatial capacity C is defined as 

K 



C = E [S*] 



E 



max , 
Re7^ A 



i=0 



(7) 



where E[-] represents the expected value. 

We can now apply properties from the point process theory to compute the spatial capacity of this class 
of Poisson networks as follows. 

Lemma 1 (average spatial capacity): The average spatial capacity C is computed for the Poisson net- 
work described in this section as 

C = \E[R% (8) 

where A is the network density and R* is the random variable that characterizes the spatial-capacity 
achievable rates over the infinite plane and network realizations. 

Proof: Let us first remind that the spatial process $ takes place in and then A ^ oo, K oo 
and R = (_Ro, -Ri, •••). Then, we proceed with the following manipulation: 



C = E 

(a) 



E 



1 

max lim — y Ri 

nan a^oo A ^ 

i=0 



A^oa A 



i=0 



XE[R* 



(9) 
(10) 
(11) 



Specifically, equality (a) considers the value of R* = {R^, R^, ...) E IZ that leads to the maximum 
spatial throughput for a given network realization, resulting in the spatial capacity S* . Equality (6) comes 
from a straight application of Campbell theorem [flT', Th. 2.2] in equation d?]), considering that $ has 
spatial intensity of A [TXs/m^]. The rates i?*'s that guarantee a successful communication are determined 
by random variables that depend on the network topology for each realization of $. As homogeneous 
PPPs are stationary and isotropic, we can apply the concept of Palm distribution and Slivnyak theorem 
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ifTTl Ch. 3] to determine the statistical proprieties of any node in $ based on a typical node. Then, we can 
determine the expected value of the spatial-capacity achievable rates over and over different realizations 
of $ based on such a node and hence the index i can be dropped, which concludes this proof. ■ 

Remark: We can generalize this results for different stationary and isotropic point processes. It is also 
possible to use a similar approach to compute the spatial capacity of more general processes using the 
Campbell theorem and Palm distributions accordingly. 

From equation ([8]), one can see that the main problem is now to derive the distribution of the spatial- 
capacity achievable rates R*, which is our focus in the next two sections. 

IV. Interference as Noise Decoding Rule 

In this section we assess the decoding rule where the receivers treat the interference as noise - or IAN 
decoders. The following corollary shows its achievable rates. 

Corollary 1 (achievable rates for IAN decoders): Assuming the noise is Gaussian and considering that 
TXs employ G-ptp codes as described in Section [III the rate Rk associated with a given link TX^-RX^ is 
achievable when IAN decoders are used if, and only if, the following inequality holds: 



Rk < log2 



Pkk 

\ jeA\{k} I 



(12) 



where A represents the set of active transmitters. 

Proof: This is a special case of ([3]) assuming that RX^ only decodes the message of TXjt while the 
other TXs are treated as noise. ■ 
We now apply this corollary to the scenario described in Section HII] to assess the spatial capacity 
of Poisson networks when receivers use IAN decoders. Before we start, however, we still need to 
characterize the propagation phenomenon. We consider here the distance-dependent path-loss modeling 
with exponent a > 2 ll29l so the channel gain between TX^ and RXj is \gij\^ = where Xij denotes 
the distance between therqj. We assume the noise power is negligible in comparison to the interference 
power (interference-limited regime). We further consider that the aggregate interference experienced by 
RXfc can be approximated by power cio related to its closest interferer. Mathematically we have the 

^This is in fact a simplified model that may lead to meaningless results for Xij < 1. As pointed in 1301 . modified versions of this model 
just increase the complexity of the analysis without providing significant differences. We can also include into our channel modeling the 
effects of random fluctuations due to shadowing and multi-path as in |12, Sec.4.1]. For our purposes, though, the incorporation of these 
phenomena only complicates the mathematical formulation without giving any further insight on the network behavior. 
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following: 1 + Yl Pkj ~ Pk,c\o- 
ie^\{fc} 

Based on these assumptions, we can derive an approximation of the probability density function (pdf) 
of the spatial-capacity-achieving rates for IAN decoders as presented next. 

Proposition 1 (approximate pdf of the spatial-capacity-achieving rates for IAN): Let R* = [Rq, ...) G 
TZ be the rate tuple that achieves the spatial capacity for the network described in Section Illl The pdf 
of i?*, V 2 G ^ over the different spatial realizations is equivalent to the pdf of a typical rate R*, 

which is approximated by 

where x > 0. 

Proof: Let us analyze a typical link TXq-RXo added to the PPP $. From Slivnyak theorem (refer to 
[fm Th. 3.1]), this inclusion does not affect the distribution of $. Without loss of generality, we assume 
that the origin of the plane is located at RXq and label the interferers TX, accordingly to their distances 
to RXq, i.e. TXi is the closest, TX2 is the second closest and so on. From our assumptions, we have 
1 + Xlfcli -Pfc ~ Pi- We then apply the path-loss model to the IAN decoder presented in equation (fT2l) . 
considering that the distances from TXq and TXi to RXq are respectively > and ri > 0, resulting in 

^0 < logs (1 + ^) , (14) 

where ri is a random variable. 

To compute the pdf of ri, we use the definition of contact zone [fTTl Defs. 1.9, 3.2] (the distance 
between a typical point and its first neighbor) to obtain the pdf of ri, resulting in [|3TI 

/,,(x) = 2A7rxe-^"^', (15) 

such that X > 0. Defining /3q = (i~°/rf" such that inequality (fT4l) still holds, then we have the following 
relation ri = c//3q° (see Fig. [T]). We now apply this variable transformation to (fTSl) and hence the pdf of 
/3q > can be obtained as 

= (16) 

° ax 

where x > 0. 

^This approximation is analyzed in f251 and it usually applied to compute lower bounds of the interference power based on dominant 
interferers 1121 . 1141 . We also discuss more about it in Section fVIll 
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Fig. 1. Illustrative example of the typical link TXq-RXo employing the IAN decoding rule, where TXi represents the closest interferer to 

* — 

RXq. To reach the highest achievable rate Rg, the relation ri — d/?,) ° must be respected such that ri is the random variable that denotes 
the distance between RXq and TXi. 



To conclude this proof, we proceed with the transformation Rq = log2(l + /Sq) remembering that PPPs 
are stationary so we can characterize any node of the network based on a typical node, dropping the index 
(refer to [Uj Sec. 3.4]). ■ 

The result just stated provides us an approximation^ of pdf of the spatial-capacity-achieving rates using 
IAN decoders over an infinite plane and over different spatial realizations of the process $. Then, we 
apply (fT3l) to approximate the spatial capacity given by ([8]), resulting in 

POO 

CiAN ~ A / xfR*{x) dx, (17) 

which does not have a closed-form solution and a numerical integration is required. For this reason, next 
we derive some proprieties|j of (flTI) that help us to understand the Cian behavior. 

Property 1 (concavity of the spatial capacity): A function /(■) is said to be quasi concave if, and only 
if, f {pxi + (1 — p)x2) > min{/(xi), /(a;2)}, where < j9 < 1. Considering that the rate that leads to 
the spatial capacity, R*, is a function of the network density A (i.e. R* = /(A)), then Cian given by (fTTl) 
is quasi concave in terms of A, where R* is a random variable characterized by the pdf ([13). 

'We discuss the tightness of the closest interferer approximation later in Section IVIII 

'^Such properties rely on the closest interferer approximation that will be discussed later on. For simplicity we hereafter refer to the 
approximate average spatial capacity as spatial capacity. 
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Proof: Let us first consider two different network densities Ai and A2 such that Ai < A2. Then, 
defining that A = pXi + (1 — p)\2 with < p < 1, we proceed with the following manipulation 

Cian(A) = (pAi + (1 - p)\2) E[/ (pAi + (1 - p)\2)] (18) 

> AiE[/(j9Ai + (l-p)A2)] (19) 
Ai E[/(Ai)] = Cian(Ai) (20) 

> A2E[/(A2)] = Cian(A2). (21) 

Inequality (a) comes from the fact that Ai < pXi + (1 —p)\2 whereas equality (6) is obtained by setting 
p = \ since the first inequality holds for all < p < 1. This proves the quasi concavity of the analyzed 
function when AiE[/(Ai)] < A2E[/(A2)]. Finally, inequality (c) is straight when AiE[/(Ai)] > A2E[/(A2)], 
which concludes this proof. ■ 

Property 2 (highest spatial capacity): The network density A* that leads to the highest spatial capacity 
given by (flTI) is obtained as the density A > which is solution to the following equation: 

j x^-^ log2(l + x)e-^"'^'^" dx = j x^'^ [X-nd^x^ - 1 j log2(l + x)e-^"'^'^'' dx. (22) 

Proof: Let us first rewrite the spatial capacity formulation using the pdf (x) given by (fT6l) . yielding 

POO 

ClAN = A / \og2{l + x)f 13* (x) dx. (23) 

^0 

Then, we recall that the Cian is quasi-concave in terms of A (Property [T]) to find its maximum value 
based on the derivative equation dCiA^/dX = 0. After some algebraic manipulation, we obtain (|22] ). which 
concludes this proof. ■ 

Property 3 (lower bound): A lower bound of the spatial capacity given by (fTTT ) is computed as 

CiAN>Aye-^-'^^(2^-i)", (24) 

where y > 0. 

Proof: To prove this property, we apply the Markov inequality as presented below: 

Pr[R* >y]< ^ E[R*] > ye-^^'^'i^'-^)" , (25) 

where Pr[i?* > y] = 1 - (x) dx and 2?^ - 1 > 0. 

Then, we multiply both sides by A, resulting in (|24l) . ■ 
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Property 4 ( upper bound): An upper bound of the spatial capacity given by (fTTT ) is computed as 

C™<Alog,(l+(3i^)"r(l + |)j. (26) 

where r(-) is the Euler gamma function defined as T{z) = t^~^e~^ dt. 

Proof: Let us apply Jensen's inequality based on the concavity of (fTTl) (refer to Property [B, yielding 

CiAN = A E[i?*] (27) 
AE[log2(l + r)] (28) 
< A log2(l + E[/3*]), (29) 

where equality (a) comes from the fact that R* = log2(l + /?*) and inequality (b) is the Jensen inequality 
for quasi-concave functions. Then, we compute the expectation of the random variable (3* using (fT6l) . 
which proves (|26] ). ■ 
Property 5 (asymptotic equivalence): Let ~ denote asymptotic equivalence of two functions, then 



CiAN-cA^-t, (30) 

TTdy ~ V ' 2/' 

Proof: To prove that two functions f{x) and g{x) are asymptotically equivalent, i.e. f(x) ~ g{x). 



when A — )• oo and c 



we should show that lim f{x)/g{x) = 1. Let us first consider the behavior of the random variable /3*, 
characterized by (fT6] ) when A — )■ oo, yielding 

lim /^.(x) =5(x), (31) 

A— >-oo 

where is the Dirac impulse function. 

This indicates that the random variable (3* tends to have the value with high probability when the 
network density increases. Now, let us consider that f3* — )■ 0, then we have the following limit 

lim !254±£) = 1 . (32) 

/3*^o /3* \n2 

Using these limits, we can manipulate the expression of the spatial capacity Cian as follows. 

lim Cian = Hm A E[log2(l + (3*)] = lim A (33) 

A— >-cxD A— J-oo A— >c» In 2 
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Fig. 2. Actual values, lower and upper bounds of the spatial capacity, Cian, versus the network density A for q = 4 and d = 1. The lower 
bound is obtained using y = 1 in equation j24b . The actual values and upper bound are computed using equations dlVt and i26\ . respectively. 



Proceeding similarly with the upper bound, we have 



E[(3*] 



lim A log2(l +E[/3*]) = lim A , ^ . 

A^oo A— >-oo In 2 

Now, we recall that the division of limits is the limit of the division, resulting in 

A E[log2(l + /?*)] 



(34) 



lim 

A^oo A log2(l + E[/3*]) 



1. 



(35) 



From this fact, we can state 



IAN 



Alog2 



(36) 



when A — )■ oo. 

To conclude this proof, we verify that 
the approximation log(l + x 



r 1 



a 



^ X valid when x << 1 into (l36l) resulting (l30l) . 



— i- when A — t- oo and then we apply 



Fig. |2] illustrates the behavior of the spatial capacity Cian and its proposed bounds as a function of 
the network density A. Firstly, one can notice that the spatial capacity has a maximum point which is 
expected from its concavity stated in Property [T] and the density A* that achieves the optimal is given 



by equation (|22l 



When densities lower than this maximum are considered, the network is spatially not 



'a closed-form solution is unknown for this equation but standard numerical methods solve it. In our case, we use Find Root from Wolfram 
Mathematica 8. 
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saturated and the spatial capacity of the network is still not in its highest value. In this situation, any 
increase of A leads to an increase of Cjan until such an inflexion point is achieved. After that point, on 
the other hand, the network spatial throughput degrades due to the proximity of the interferers, strongly 
reducing the average of the link rates R*. Consequently, Cjan becomes a decreasing function of A. 

From Fig. |2l we can also evaluate the proposed upper and lower bounds of the spatial capacity. As 
one can notice the lower bound proposed in Property |3] is loose, regardless of A. In fact, the main use 
for this bound is to prove the relation between the spatial capacity and the maximum spatial throughput 
achieved with fixed rates, as discussed later on. Regarding the upper bound introduced in Property S] 
when A increases, the upper bound become tighter and tighter, as predict by Property |5] In other words, 
the upper bound has the same value as the spatial capacity Cian when A — )• oo as shown in Fig. |2l 

In the next section, we apply the same approach used here to derive the approximate spatial capacity 
and its properties when OPT decoders are considered. 



As previously discussed, the optimal decoding strategy when Gaussian point-to-point codes are used in 
wireless networks with multiple transmitter-receiver pairs consists in jointly decoding some few messages 
from the strongest interferers, while the rest is treated as noise. Based on this result, we obtain the 
achievable rates for links whose receivers use OPT decoders as follows. 

Corollary 2 (achievable rates for OPT decoding rule): Assuming that the noise is Gaussian and con- 
sidering that TXs use Gaussian point-to-point codes as described in Section [III then the rate Rk associated 
with a given link TX^-RX^ is said to be achievable when the OPT decoder is employed if, and only if, 
the following inequality holds: 



where Al represents the subset of transmitters whose messages are decoded by receiver k and AlUAl = A 
is the set of all active transmitters in the network. 

Proof: To obtain (l37l) we proceed with a simple manipulation of equation ([3]), isolating the rate Rk 



As in Section |IVl we apply the theorem stated above to statistically characterize the achievable rates 
over different spatial realizations using now the OPT decoding rule and then assess the average spatial 
capacity of the network described in Section Hill which is given by equation ([7]). For the scenario used 



V. Optimal Decoding Rule 




(37) 



related to TX^-RX^ link by considering the subsets Al that lead to achievable rates. 
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throughout this section, however, the analysis is more complex since the receiver node should choose 
the subset of messages that will be jointly decoded and then verify whether the coding rate of its own 
transmitter is achievable, given all other coding rates. By construction, all receivers proceed in the same 
way and hence the analysis becomes a very intricate combinatorial problem. For this reason, in order to 
derive the pdf of the achievable rates for the OPT decoders, we resort to some approximations that will 
be justified afterwards. 

As before, we only consider the deterministic path-loss (refer to Section HVT) and that the sum of the 
interfering signals observed by RX^ that are treated as noise can be approximated by the signal from the 
closest interferer amongst them, whose power is denoted Pfc.cio- If the noise power is negligible compared 
to Pfc,cio, then 1 + Yl Pkj ~ -Pfc.cio- Based on these assumptions, we can state the following proposition. 

Proposition 2 (approximate pdf of the spatial-capacity-achievable rates for OPT): Let us denote the rate 
tuple that achieves the spatial capacity for the network described in Section UlI] as R* = (_Rq, Rl, ...) G TZ. 
Then, the pdf of Rl,'i k E A follows the pdf of a typical rate R*, denoted by fR*{x) and approximated 

^ J^^^^rrf e-^'^'^Vii>(x|n = t) (38) 



TU) 

where fR*\n{x\n) is the pdf of R* given that 1 + n messages are jointly decoded and is approximated by 
/«.|„(x|n)^ln4 [-TT^) ' ^ ^' ^^^^ 

such that X > 

l+n 

Proof: Let us first deal with the typical link TXq-RXq. Without loss of generality, we place the 
origin of the Cartesian plane at RXq and assume that all nodes that are closer to RXq than TXq have 
their messages jointly decoded with TXq message (see Fig. [3]). From the distance-dependent path-loss 
modeling, the closer the TX, the higher the power, and then this choice of the subset is justified. 

For each network spatial realization, we consider that a number n associated with the transmitters whose 
messages are decoded by RXq is known, which yields the following inequality 

V -T O,clo / V -T O,clo / 



One can observe from (|37T ) and (|40] ) that rate tuples that satisfy Rq + J2^=i < log (^^ + 
are always achievable. 



(l+n)Poo 

fo.clo 
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TXi 



Fig. 3. Illustrative example of the typical link TXq-RXo employing the OPT decoding rule. The blue TX has its message jointly decoded 
with TXo message and TXi is the closest interferer to RXo whose signal is treated as noise. The random variable ri denotes the distance 
between RXo and TXi such that ri > d. 



Defining /3q = Poo/-Po,c1oj we use similar steps to the ones used in the proof of Proposition |2l but 
considering now that ri > to compute the pdf fp* (x) as 



where x > 1 and fi3*{x) = when x < 1. 

Then, we assume that Rq + X;r=i -Ri ~ (1 + n)RQ to obtain (1 + n)RQ = log (1 + (1 + n)l3^). By 
applying such a transformation, we can find the pdf of Rq given n. Here we use the proprieties of Palm 
distributions in PPPs (refer to the proof of Proposition |2j and then the index can be dropped, resulting in 



Remark: In addition to the closest interferer treated as noise approximation, this proposition is based on 
other two main assumptions: (i) the detected power at RXq related to the 1 + n jointly decoded messages 
is equal to (1 + n)Poo and (ii) the sum rate associated with those messages is given by (1 + n)Ro. In 
assumption (i), we use the lower bound given by (|40] ). which indicates that we underestimate the aggregate 
power, while in (ii) we approximate the sum of 1 + n random variables that follows the same distribution 
by one random variable multiplied by 1 + n. We argue that the underestimation byproduct of assumption 
(i) leaves us some room for variations in the sum rate approximation used in (ii). In addition, due to 
the homogeneity of the spatial process, Rq + XlILi -Rj ~ (1 + n)Ro leads to a reasonable approximation. 
Simulations results are presented in Section IVIII where we discuss such approximations. 




(41) 



equation (|39l) . To unconditioned the pdf /i?*|„(a;|n), we apply the definition of Poisson process to compute 
the probability that n = i points lie in the area 7r(f, concluding this proof. ■ 
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Here we can compute the approximate average spatial capacit>o Copt when the OPT decoding rule is 
employed as 

POO 

Copt- A / xfu^ix) dx, (42) 



where fR*{x) is given in Proposition [2l 

The integral in (l42l) is analytically unsolvable (we can rely on numerical solutions, though). To gain 
more insights on the system performance, we next derive some properties of the spatial capacity. 

Property 6 (concavity): Considering that the rate R* is a function of the network density A, then Copt 
given by (|42l) is quasi concave in terms of A, where R* is a random variable given by (|38] ). 

Property 7 (lower bound): A lower bound of the spatial capacity given by (|42] ) is computed as 



2_ 

^2\i x„.2/'2(l+')!'-r 



Copt > A > ye V / , (43) 

^ r z 



where v > '"^2(^+0 foj- all i > 0. 

^ i+i — 

Property 8 ( upper bound): A upper bound of the spatial capacity given by (|42] ) is computed as 

S ^ E ^ ^ '08^ (1 + (1 + * r (1 + ^. A..=) e-) . (44, 



where r(-, ■) is the incomplete Gamma function, which is defined as T(z,a) = * dt. 

Property 9 (asymptotic equivalence): Let ~ denote asymptotic equivalence of two functions, then 



Cop. - A I £^ 10., (1 . (1 . (3^) ' r (1 + ^, A.^') e^^) , (45, 

when A — )■ oo. 

The proofs of these properties follow the same principles used in the previous section so we do not 
present them here. It is worth pointing though out that the proofs of (|43])-(|45]) begin by assuming that 
the number 1 + n of jointly decoded messages is known. Then, we use the fact that the unconditioned 
spatial capacity is a linear combination of the conditioned capacities with weights given by the Poisson 



probabilities that n = i nodes lie in a area of nd"^, resulting in the probability ^ e 



Next, we present a proposition that states that the OPT decoding rule always leads to a better perfor- 
mance than the IAN rule. 

'"As in the previous section we use tlie term spatial capacity to refer to the approximate average spatial capacity. 
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Proposition 3 (Cqw vs. CianJ-' For any given network density A, 

Copt > Cian- (46) 

Proof: To prove this proposition, we first use the identity alog2(l + x) < log2(l + ax) V a > 1 to 
verify that 

r.i- + + > r.i-^ log,(l + .) (47) 

Ji 1 + n Jq 

where n G N represents the number of jointly decoded messages by the typical link XXq-RXq. 

Recalling that all TXs that are located closer to RXq than TXq have their messages jointly decoded by 
RXq, then the number n follows a Poisson distribution over an area of 'Kd?. Using this fact and multiplying 
both sides of (l47l) by we have the following inequality: 

V i^e--^ r forix) + + > V (^e--^ r log.d + X) dx, 

T{{) A 1 + ^ "t^ r(z) Jo ^ 

(48) 

where f^F^ix) and /^^^^{x) are given by (1411 and (fT6l) . respectively. 

To conclude this proof, we multiply both sides of (|48T ) by A, obtaining then (|46|) . ■ 

In Fig. m we present the spatial capacity Cqpt given by (l42l) as a function of A together with its proposed 
upper and lower bounds. One can observe that the lower bound given by Property |7] is very loose for the 
value of the constant y that was arbitrarily chosen {y = 2). This bound, however, can be improved by 
tuning the constant y in accordance to the number of jointly decoded messages. Such an improvement 
in the proposed bound will be discussed in the next section when we apply it to analytically assess the 
performance of networks where predetermined fixed rates are imposed. 

Turning your attention to the values of Cqpt given by (|42l) . one can easily see that it is an increasing 
function of A. For lower densities, Cqpt increases faster since the probability that an interfering TX has 
its message jointly decoded is also low and, consequently, the rate is constrained by the interferers that 
are treated as noise, indicating that Cqpt is limited by the low spatial reuse. When A increases, on the 
other hand, more messages from interfering TXs start being jointly decoded, which diminishes the Cqpt 
rate of increase. Furthermore, we can observe that the upper bound proposed in Property [8] is a good 
approximation to Cqpt for all densities A especially when A — )► oo, corroborating Property ^ 

By comparing the results shown in Fig. |2](IAN) and Fig. |4](0PT), one can see that the OPT decoding 
rule provides higher spatial capacities, regardless of the network density, as predicted by Proposition |3] 
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The performance gain obtained with the OPT decoder in comparison with the IAN decoder indicates that 
the mechanism of joint detection used here is a good way to deal with the strongest interferers. A more 
detailed comparative analysis between OPT and IAN decoding rules is presented later 

In the following section, we compare the results obtained so far with the most usual approach found 
in the literature (e.g. [fT2l '): coding rates are fixed for a given network density and set to optimize the 
average spatial throughput regardless of a specific network topology, leading to outage events (i.e. some 
pairs use coding rates above their channel capacity). 

VI. Spatial Throughput Optimization Using Predetermined Fixed Rates 

We now focus our attention on scenarios where TXs set their coding rates to the fixed values that 
leads to the highest average spatial throughput, given that the TXs are aware of how many messages are 
jointly decoded by their RXs. Using this scheme, groups of TXs use the same fixed coding rates and 
then an optimization problem is formulated to find these rates such that the average spatial throughput 
is maximized. As a consequence the optimal choice of coding rates, as discussed later on, is not in the 
network capacity region, stated in Theorem [T] leading to outage events for some links. Next, we present 
the definition of the aforementioned optimization problem. 

Definition 5 (highest spatial throughput): The spatial throughput optimization problem for a network 
where TXs have fixed coding rates is defined as 

r = maxE[5], (49) 
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where T is the maximum spatial throughput, R = (i?o, -Ri, •••) represents the set of fixed coding rates Ri 
used by the TXs such that i is the number of jointly decoded messages in addition to the desired one, and 
S is the spatial throughput given by Q, where only the successful transmissions are taking into account. 

When the IAN decoding rule is used, there is no jointly decoded message and then the optimization 
is only related to one fixed coding rateLlJ. We now present two propositions that state the highest spatial 



reQ. 



throughputs for IAN and OPT decoders applying the network modeling used beforq 

Proposition 4 (highest spatial throughput for IAN): The highest spatial throughput Tian achieved when 
IAN decoders are used is given by 

TIan = Alog2(l + /3*)e-"-'^'^*", (50) 

where /3* is the value of /3 > 0, which is solution of 

P= (^^X7id\l + /3)ln{l + (3)y'\ (51) 

Proof: Let us first write the spatial throughput given by dSj) for this scenario as 

S = XRP„ (52) 

where R is the fixed coding rate used by all TXs and is the success probability associated with R. 

We proceed here similarly to the proof of Proposition [Hand then apply the relation R = log2(l + /?), 
where R, (3 > 0. From the dominant interferer assumption, an outage event occurs whenever an interfering 
TX node lies inside the area defined by the circumference centered at the RX node and with radius d(3^. 
Using the Poisson distribution, we have that Pj = e"^'^^'^^" . Hence, we can rewrite equation (|52] ) as 

5 = Alog2(l + /3)e-^"'^'^", (53) 

which is a concave function of f3. 

Hence, we compute (3* which is the solution of the derivative equation dS /d(3 = 0, resulting after some 
manipulation in (|5TI) . To conclude this proof, we use /3* into equation (|53l) . obtaining (|50l) . ■ 

Proposition 5 (highest spatial throughput for OPT): The highest spatial throughput Topt achieved when 
OPT decoders are used is given by 

r„„ ^ A f ^ ,og,(l + (1 + m) (54) 

where, for each i = 0, 1, 2, /3* is found as the value of /3i > 1, which is solution of 

"This is the usual approach as in 1121 . 11241 . 

'^Once again we use the closest intereferer treated as noise approximation. Besides, the term highest spatial throughput refers to the highest 
approximate spatial throughput. 
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V(l + z)a J 

Outline of proof: To prove this proposition, we follow the same steps used in the proof of Proposition 
in considering these basic differences: /3j = (i^"/rf" > 1 (since messages from TXs closer to a given 
RX than its own TX are jointly decoded and then ri > d) and the optimization is proceeded for each 
i = 0,1,2,... which yields (1551) . To conclude this outline, we average the spatial throughputs by the 
Poisson probabilities that i nodes lie in the area ncP, resulting in (l54l) . ■ 

Here we apply Properties |3] and |7] to obtain an analytical relation between the spatial capacity C and 
the highest spatial throughput T using fixed rates for either decoding rules. 

Proposition 6 (C vs. T): For a given density A and assuming that all links use the same decoding rule 
(either IAN or OPT), then 

C>r. (56) 
Proof: This statement is a consequence of Property |3l when we set the constant y = log(l + /?*) in 
(l24l) . yielding (l50l) . Similarly, we use Property Ul applying for each different z G N a different constant y 
in (1431) such that = i°s(^+(|+')^i ) ^ which yields (l54l) . concluding this proof. ■ 
Fig. [5] shows the maximum spatial throughput following the formulation derived in this section. As 
proved in Proposition |6l T is always lower than or equal to C for the same density and the same decoding 
rule. This is justified by the methodology used to derive the spatial capacity, which allows for a coding 
rate setting based on the relative positions of the nodes for each different realization of the spatial process. 
When fixed rates are used, the transmitters code their messages using a fixed rate that depends only on 
the number of other messages that are jointly decoded by their own receivers. By optimizing based only 
on the average behavior of the network, some RXs cannot successfully decode their messages for specific 
topologies, which decreases the spatial throughput. Therefore, the spatial-capacity-achieving strategy has 
always a better performance. Besides given the decoding rule employed, the curves of T and C have a 
similar shape. 

Fig. [5] also shows that the spatial capacity obtained when OPT is used has a huge gain if compared 
with the IAN option. As predicted by Proposition |3l this result reflects that the OPT rule is able to avoid 
the strongest interferers by jointly decoding their messages. When the density A is low, both OPT and 
IAN decoders have approximately the same performance since the probability that a interferer is closer 
to a given RX than its own TX is very low. Increasing A, such a probability also increases and then the 
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Fig. 5. The highest spatial throughputs T using fixed coding rates given by dSOt and i54\ . and the spatial capacities C given by l ll7t and 
l l42b as a function of the network density A for IAN and OPT decoding rules, d — 1 and a — 4. 

differences between the strategies become apparent as the closest interferer is the limiting factor for IAN, 
while such node may have its message jointly decoded when OPT is used, what decreases the harmful 
effects of the nearby interferers. 

VII. Discussions 

So far we have showed that, for same network density, (i) OPT decoders outperform IAN, and (ii) the 
spatial-capacity-achieving strategy outperforms the average spatial throughput optimization when receivers 
employ the same decoding rule. Nevertheless we still need to discuss some possible limitations of our 
finds, namely the tightness of our approximations and the feasibility of each decoding rule for practical 
implementations. In the following subsections we deal with both aspects, identifying why our results 
are important even when our approximation is poor and for which circumstances the design setting that 
provides the worst performance is more suitable than the optimal. 



A. Tightness of our approximation 

Here we discuss the validity of the "closest interferer treated as noise approximation" used to derive the 
approximate performance of both decoding rules. Figure |6] shows the average spatial capacity C computed 
using our analytical approximation and Monte Carlo simulation as a function of the network density A 
for both decoders!^. For both IAN and OPT, the lower the density is, the better our approximation works. 

'^The results for the highest average spatial throughput presented in Section IVll follows the same trends and thus we exclude them from 
the paper. 
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Conversely, increasing the density, our approximate spatial capacity gets looser and looser. 

The closest interferer approximation is in fact a lower bound of the aggregate interference [[T2l . leading 



then to an upper bound of the actual spatial capacity .^^ 
equivalent to the actual values when A — > [fT2ll . ES 



his bound have been proved to be asymptotically 
For higher densities, the closest interferer treated 



as noise tends to contribute less to the aggregate interference experienced by the receivers, worsening our 
approximation. Besides, we obtained our numerical results using the path-loss exponent a = 4 and Weber 
et al. showed that lower exponents lead to looser bounds 09]. 

Yet, the comparison between the IAN and OPT decoders is fair since the results presented in Sections 
UV] and |V] rely on the same approximatioi|lf|. We further argue that our approximation has no effect in the 
trade-off analysis done in this paper and Figure |6] illustrates this fact by showing that the OPT always 
outperforms IAN in similar scales: the ratios Cian/Cqpt obtained via simulation or via our approximations 
have similar values when considering the same A. As the proposed formulation provides a much simpler 
way to assess the network performance than numerical simulations, we reinforce the contribution of this 
paper even when our approximations provide poor bounds. 

All in all, we believe that our main messages - OPT is better than IAN, and spatial-capacity achieving 
strategy is better than the best fixed rate one - are unaffected by our approximations. Despite of these facts, 

''*In our point of view this asymptotic analysis is unsuitable for the study carried out here; we assume an interference-limited network, 
which opposes the idea of very low density of interferers. When A — >■ 0, we see the network in its noise-limited regime. 
'^We can argue in the same way to say that the analysis presented in Section IVTl is also fair. 
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the optimal strategy is many times unfeasible for practical implementation as discussed in the following. 



B. Design setting and mobility pattern 

Throughout this paper we have shown that the best design option in terms of spatial throughput is 
to employ OPT decoders and apply the spatial-capacity achieving scheme. This solution, however, has 
drawbacks: (i) RXs require the knowledge of the codebooks of the jointly decoded messages and (ii) OPT 
decoders are computationally more complex than IAN. 

Knowing them, we argue that the use of either/both OPT and/or spatial-capacity achieving strategy 
is unfeasible for (highly) mobile topologies. Under this topology, the neighbors of any given receiver 
vary very fast, making impossible the joint decoding procedure. Shopping malls and streets where people 
frequently come and go can exemplify this scenario. If this is the case, even though the design setting 
employing IAN decoders with fixed rate optimization is far from the optimal performance, it is a more 
suitable choice. 

Conversely, when (quasi-)static networks are considered, the optimal strategy becomes possible. In this 
case, receiver nodes must known the codebooks of their strongest interfering nodes and jointly decode their 

messages. In addition, the links must coordinate their coding rates to be in the network capacity region. 
Smart homes, industry plants and other kind of machine-to-machine communications can exemplify this 
mobility pattern. 

Besides, there are other aspects that may be prohibitive for OPT. For instance, many applications 
require secrecy and then the codebook knowledge makes OPT unfeasible even for static topologies. Other 
applications need fast processing time, which is also unfeasible when many interfering messages are 
jointly decoded. Anyway, this dependence of the topology must be taken into account when the network 
is designed. Furthermore, the mobility pattern of the network can also change over time - for example, 
offices during the night are quasi-static, while highly dynamic during parts of the working hours. All 
these aspects indicates the needs for adaptive algorithms that may assess the network state and proceed 
with their optimization accordingly. 

VIII. Final Remarks 

In this paper we studied the spatial throughput of interference networks using Gaussian point-to-point 

codes. More specifically we evaluated two different decoding rules: (i) treating all interfering messages as 
noise - IAN, and jointly decoding the messages whose detected power is higher than the desired message 
power while treating the others as noise - OPT. 
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We proposed an approximation of the highest spatial throughput for Poisson distributed networks such 
that all links work in their capacity region (spatial-capacity achieving strategy). We then stated several 
properties of our approximation using either decoders and prove that, when the same network density is 
assumed, (i) the OPT always outperforms IAN, and (ii) the spatial-capacity-achieving strategy is always 
better than the predetermined fixed rate optimization, where transmitters code their messages in order to 
optimize the average spatial throughput. 
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